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ABSTRACT 

We  present  a  new  method  and  performance  data  related  to  volumet¬ 
ric  acoustic  intensity  imaging  using  a  spherical  microphone  array 
(SMA)  consisting  of  a  dual,  concentric  rigid  and  open  SMA.  The 
dual,  concentric  array  was  designed  to  improve  the  frequency  range 
of  a  standard  SMA.  We  apply  standard  techniques  associated  with 
interior  spherical  near-field  acoustic  holography  (NAH)  and,  in  par¬ 
ticular,  consider  issues  related  to  the  optimal  use  of  information  from 
both  arrays  for  NAH  projection  and  the  advantages  that  thus  accrue 
from  utilising  a  dual,  concentric  SMA. 

Index  Terms —  Array  signal  processing,  Microphone  arrays, 
Sound  field  recording,  Spherical  microphone  array,  Near-field 
acoustic  holography 

1.  INTRODUCTION 

Microphone  array  techniques  associated  with  near-field  acoustic 
holography  (for  review  see  [1])  and  discrete  inverse  acoustic  prob¬ 
lems  (for  review  see  [2, 3])  have  been  well  developed  over  the  past 
several  decades.  Recently,  there  has  been  considerable  development 
of  interior  spherical  NAH  using  SMAs  (see,  e.g.,  [4-7]).  SMAs 
present  a  wondrous  opportunity  to  characterise,  analyse,  and  recre¬ 
ate  three-dimensional  sound  fields,  which  has  probably  led  to  the 
wealth  of  recent  research.  We  have  previously  presented  a  dual, 
concentric  rigid  and  open  SMA  comprised  of  64  microphones,  that 
is  still  compact,  and  which  increases  the  available  frequency  range 
of  analysis  (see  [8]).  In  this  paper,  we  examine  the  consequences  of 
employing  a  dual,  concentric  SMA  for  volumetric  acoustic  intensity 
imaging.  That  is  to  say,  the  holographic  projection  of  the  acoustic 
intensity  field  surrounding  the  SMA  can  be  computed  using  a  com¬ 
bination  of  interior  NAH  analysis  applied  to  both  the  open  and  rigid 
arrays  and  we  explore  a  method  that  combines  information  from 
both  arrays  in  an  optimal  manner.  We  have  not  previously  seen  in 
the  literature  the  application  of  this  method  of  NAH  analysis  to  a 
dual,  concentric  rigid  and  open  SMA.  We  present  real  measurements 
and  field  data  that  demonstrate  the  performance  capababilities. 

2.  METHODS 

A  dual,  concentric  rigid  and  open  SMA  was  constructed  by  placing 
a  rigid  SMA  at  the  center  of  a  larger  open  SMA  (see  Figure  1). 
The  rigid  SMA  was  machined  from  a  Nylon  plastic,  Ertalon,  into 
two  hemispherical  halves  resulting  in  a  sphere  with  radius  1.63  cm. 
The  open  array  was  constructed  using  horizontal,  circular  steel  rings 
mounted  on  a  vertical  steel  ring,  resulting  in  a  sphere  of  radius 
6.0  cm.  Both  spherical  arrays  were  attached  to  and  supported  by  a 
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vertical,  metallic  rod.  Thirty-two  DPA  type  4060-BM  omnidirec¬ 
tional  microphones  were  arranged  on  both  the  rigid  and  open  spheres 
using  a  spherical  7-design  optimised  for  discrete,  spherical  integra¬ 
tion  (described  in  [9]).  The  relative  orientation  of  the  rigid  and 
open  SMAs  was  left  flexible.  The  microphone  cables  were  neatly 
run  along  the  support  rod  for  the  inner  SMA  and  along  the  circular 
steel  rings  for  the  outer  SMA.  While  it  seems  intuitively  reasonable 
to  distribute  half  of  the  sixty-four  microphones  to  both  SMAs  in 
order  to  maximise  the  available  frequency  range,  we  have  shown 
mathematically  that  this  actually  is  the  case  (for  details  see  [8]).  In 
particular,  we  designed  the  dual,  concentric  SMA  so  that  is  has  a 
minimum  threshold  signal-to-noise  ratio  of  30  dB  when  used  as  a 
constant  third-order  spherical  harmonic  beamformer  over  a  900  Hz 
to  16.0  kHz  frequency  range.  The  directivity  of  the  system  remains 
approximately  constant  across  the  frequency  range  by  analyzing 
frequencies  between  900  Hz  and  2.84  kHz  using  the  open  SMA 
and  frequencies  between  2.84  kHz  and  16  kHz  with  the  rigid  SMA 
(see  [10]). 

The  performance  of  the  dual,  concentric  SMA  was  analyzed  by 
simultaneously  measuring  the  acoustic  tranfer  functions  for  all  64 
microphones  for  each  of  393  directions  around  a  sphere  of  radius 
1  meter.  The  measurements  were  performed  in  an  anechoic  chamber 
with  a  loudspeaker  mounted  on  a  robotic  arm  that  could  rotate  about 
the  center  of  the  chamber  with  a  precision  better  than  one  degree 
of  spherical  angle.  The  acoustic  transfer  functions  were  measured 
using  a  log  sine  sweep  (8  s,  100  Hz  to  20  kHz).  To  make  the  record¬ 
ings  we  used  eight  8-channel  Pro  Tools  Digidesign  Pre  preampli¬ 
fiers  with  digital  gain  settings  and  four  16-channel  Apogee  AD-16X 
analogue-to-digital  converters  (48  kHz  sample  rate,  24-bit  resolu¬ 
tion).  The  ADAT  output  from  the  analogue-to-digital  converters  is 
sent  to  a  RME  ADI-648  MADI-to-ADAT  converter,  which  converts 
the  eight  8-channel  ADAT  signals  to  one  64-channel  MADI  signal. 
The  64  microphone  signals,  in  MADI  format,  are  recorded  onto  a 
personal  computer  (PC)  using  a  RME  Hammerfall  DSP  MADI  PCI 
sound  card.  The  PC  is  a  standard  workstation  computer  with  an  Intel 
Core  2  Duo  CPU  running  at  2.4  GHz  and  with  2  GB  of  RAM. 

2.1.  NAH  using  a  dual,  concentric  SMA 

The  primary  issue  addressed  in  this  paper  is  how  to  use  a  dual  SMA 
appropriately  for  volumetric  imaging  using  NAH.  In  other  words, 
with  the  dual,  concentric  SMA  there  are  two  measurement  surfaces 
and  one  would  like  to  combine  information  optimally  and  robustly. 
As  the  spherical  coordinate  system  is  separable,  NAH  proceeds  by 
first  determining  the  holographic  projection  coefficients  based  on  the 
Fourier-Bessel  representation.  The  projection  coefficients,  Pmn(k), 
are  determined  from  the  acoustic  pressure  measurements,  p(kr ,  Qj) 
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Fig.  1.  A  photo  is  showing  the  dual,  concentric  SMA  located  inside 
the  anechoic  chamber  in  which  its’  performance  was  measured.  The 
moveable  loudspeaker  is  shown  to  the  left  of  the  SMA. 


as  shown  below: 
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where,  k  =  ^ ,  is  the  wave  number,  A,  is  the  wavelength,  r,  is  the  ra¬ 
dius  at  which  the  sound  pressure  is  measured,  a  <  r,  is  the  radius  of 
the  inner  rigid  spherical  baffler,  Qj  =  ( Oj ,  4>j),  is  the  angular  posi¬ 
tion  of  microphone  j  on  the  spherical  microphone  array  in  spherical 
coordinates,  M  =  32,  is  the  number  of  microphones  on  either  the 
open  or  rigid  SMA,  aj  —  jj  ,  is  a  weight  for  each  position  related 
to  the  spatial  sampling  scheme,  p(kr,  Qj),  is  the  pressure  recorded 
by  the  microphone  at  position  Qj  on  the  spherical  microphone  ar¬ 
ray,  bn(kr ,  ka)  is  a  modal  coefficient,  and  Y™(-)  are  the  spherical 
harmonic  functions.  The  Y™(-)  are  defined  as 


VT(«)  =  rnm(M)  =  y/2n4^1^  +  ^|jT(cosg)eim^  (2) 

where,  P™(-),  is  the  associated  Legendre  function.  Note  that  com¬ 
plex  conjugation  is  denoted  by  an  asterisk.  For  the  dual,  concentric 
SMA,  bn(kr ,  ka),  is  defined  as 

bn(kr,  ka)  =  jn(kr)  -  hn(kr),  (3) 

where  k,  r  and  a  are  as  defined  previously,  jn(- ),  j'n(')>  hn(-), 
and  h'n(‘)  are  the  spherical  Bessel  function,  the  spherical  Hankel 
function  of  the  second-kind,  and  their  derivatives,  respectively,  and 
i  =  y/—l. 

Once  the  holographic  projection  coefficients  are  determined  one 
can  project  or  propagate  the  sound  field  to  a  new  radius  using  the 
following  equations: 

N  n 

P(kr,Q o)  =  J2P^(k)jn(kr)  £  Y™(Q 0),  (4) 
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where  P(kr ,  Oo)  is  the  pressure,  vr(kf ,  O0)  is  the  radial  velocity, 
ve  ( kr ,  Oo)  is  the  velocity  in  the  0-direction,  v $  ( kr ,  Oo)  is  the  veloc¬ 
ity  in  the  ^-direction,  r  is  the  projection  radius,  Oo  is  the  projection 
angle  in  spherical  coordinates,  p  is  the  density  of  air,  c  is  the  speed 
of  sound,  N  is  the  truncation  order  for  the  Fourier-Bessel  expansion 
and  the  other  terms  are  as  defined  previously.  Note  that  the  above 
equations  differ  slightly  from  those  given  in  [1]  because  we  have  ex¬ 
pressed  ve(kr ,  fio)  and  v^i^kr,  fio)  in  radians.  The  acoustic  inten¬ 
sity  field  is  given  by  (Pv*).  The  equation  given  for  the  calcu¬ 
lation  of  the  holographic  projection  coefficients,  Pmn(k),  highlights 
the  fact  that  the  pressure  measurements  are  scaled  by  the  modal  co¬ 
efficients,  bn(kr ,  ka).  This  scaling  removes  the  dependence  of  the 
projection  coefficients  on  the  radius  of  the  measurement  surface.  In 
other  words,  the  Pmn(k)  should  in  some  sense  be  independent  of 
the  measurement  surface  and,  in  particular,  independent  of  whether 
the  rigid  or  open  measurement  surface  is  used  for  their  calculation. 
However,  this  is  not  actually  the  case  with  real  measurements  be¬ 
cause  the  modal  coefficients  are  different  for  the  two  measurement 
surfaces  (see  Figure  2).  The  physical  implications  of  this  difference 
are  that  for  a  given  sound  field  the  modal  coefficients  provide  a  mea¬ 
sure  of  the  signal  level  that  is  present  at  the  microphones  on  the 
measurement  surface  for  a  given  order  of  the  Fourier-Bessel  repre¬ 
sentation.  Thus,  we  propose  a  simple  and  effective  means  for  cal¬ 
culating  the  holographic  projection  coefficients:  Use  the  projection 
coefficients  which  are  associated  with  the  measurement  surface  with 
the  larger  modal  coefficient  for  a  given  order  of  the  Fourier-Bessel 
expansion. 


2.2.  Reguiarisation,  Spatial  Aliasing,  and  Measurement  Noise 

The  impact  of  spatial  aliasing  error,  microphone  placement  error  and 
measurement  noise  error  on  the  performance  of  the  dual,  concentric 
SMA  is  similar  for  both  beamforming  and  NAH.  In  designing  the 
dual,  concentric  SMA  we  made  extensive  theoretical  calculations  of 
these  errors  (see  [8])  to  ensure  that  the  array  would  operate  at  third- 
order  for  the  specified  frequency  range  (900  Hz  to  16.0  kHz).  The 
description  of  array  errors  differs  slightly  for  beamforming  [11]  and 
NAH  [7],  but  the  mathematical  essence  is  the  same.  Spatial  aliasing 
error  limits  performance  at  the  higher  frequencies,  while  measure¬ 
ment  noise  error  limits  performance  at  the  lower  frequencies.  Spa¬ 
tial  aliasing  error  arises  when  the  Fourier-Bessel  expansion  used  for 
the  holographic  projection  (refer  to  Eqn.’s  4-7)  is  truncated  to  a  fi¬ 
nite  order  (N).  This  truncation  provides  a  method  to  regularise  the 
holographic  projection  solution  and  is  similar  to  the  k- space  filter 
commonly  used  in  planar  NAH  (see  [1]  for  details).  In  [1 1],  Rafaely 
derives  the  spatial  aliasing  error  and  the  measurement  noise  error 
for  a  SMA  in  a  plane-wave  soundfield  with  the  array  looking  in  the 
direction  of  the  plane- wave,  as  is  common  in  the  literature.  Using 
Rafaely’ s  error  formulation  we  demonstrate  in  [10]  that  our  dual, 
concentric  SMA  is  accurately  described  by  these  errors  in  practice 
as  well  as  in  theory. 
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Fig.  2.  The  amplitude  of  the  modal  coefficient  for  the  inner  rigid 
SMA,  and  the  outer  open  SMA  is  shown  for,  (a),  orders  zero  and 
one  and,  (b),  orders  two  and  three. 


3.  RESULTS 

In  Figure  3  we  show  the  acoustic  intensity  field  that  is  obtained  when 
NAH  is  naively  applied  to  the  dual,  concentric  SMA  for  a  2600  Hz 
sound  source  located  on  the  x-axis  at  1  m  from  the  center  of  the  ar¬ 
ray  (as  described  in  Section  2).  By  naively,  we  mean  that  we  do 
not  combine  projection  coefficients  from  the  two  SMAs.  We  show 
a  holographic  projection  of  the  acoustic  intensity  field  as  projected 
outward  from  the  inner  rigid  array  to  a  2  cm  projection  radius.  We 
also  show  a  holographic  projection  of  the  acoustic  intensity  field  as 
projected  inward  from  the  outer  array  to  a  5  cm  projection  radius. 
Note  that  all  holographic  projections  described  in  this  section  are 
truncated  to  third  order.  For  this  first  example,  we  have  purpose¬ 
fully  selected  a  frequency  at  the  boundary  of  which  either  the  rigid 
or  open  array  can  operate.  What  we  would  like  to  highlight  is  that 
even  at  these  relatively  small  projection  distances,  we  do  not  obtain 
a  correct  holographic  projection.  For  example,  the  field  projection 
using  the  open  array  clearly  has  too  much  curvature.  We  can  start 
to  correct  and  compensate  for  these  defaults  by  selecting  a  combina¬ 
tion  of  projection  coefficients  from  the  inner  and  outer  arrays  based 
on  the  modal  coefficients  as  described  in  Section  2.1.  So,  for  ex¬ 
ample,  consider  now  Figures  4(a)  and  4(b)  which  constrasts  naive 
projection  with  “combination  projection”.  In  these  figures,  we  show 
holographic  projection  of  the  acoustic  intensity  field  at  a  projection 
radius  of  3.13  cm  (the  geometric  mean  of  the  inner  and  outer  array) 
for  a  2850  Hz  sound  source  located  on  the  x-axis  at  1  m.  For  the  com¬ 
bination  projection,  the  order-0  and  order- 1  projection  coefficients 
are  determined  from  the  rigid  array,  while  the  order-2  and  order-3 
coefficients  are  determined  from  the  open  array  (refer  to  Figure  2). 
Each  figure  actually  shows  two  projections,  one  for  the  combination 


Fig.  3.  Naive  holographic  projection  is  shown;  details  as  in  text. 

projection  and  one  for  the  naive  projection.  Clearly,  the  combination 
projection  gives  better  performance  and  allows  the  array  to  oper¬ 
ate  in  a  region  it  would  not  normally  be  able  to  do  so.  Similarly, 
Figures  5(a)  and  5(b)  show  a  holographic  projection  at  a  projection 
radius  of  2.0  cm  for  a  5000  Hz  source  located  on  the  x-axis  at  1  m. 
At  this  frequency  for  the  combination  projection,  the  order-0,  order- 
1,  and  order-2  coefficients  are  determined  from  the  rigid  array  and 
the  order- 3  coefficients  are  determined  from  the  open  array  (refer  to 
Figure  2).  Again,  we  see  a  clear  improvement  in  the  quality  of  the 
holographic  projection. 

In  order  to  quantify  the  above  results,  we  used  the  rigid  and  open 
arrays  to  estimate  the  free-field  sound  pressure  field  (i.e.,  the  pres¬ 
sure  field  without  the  array  present  and  with  no  scattering  from  the 
rigid  spherical  baffle)  at  their  corresponding  measurement  surfaces 
at  the  32  microphone  positions  and  treat  this  as  ground  truth.  We 
then  used  the  projection  coefficients  from  the  rigid  array  alone  and 
the  open  array  alone  to  estimate  the  pressure  field  at  both  measure¬ 
ment  surfaces  and  compare  this  with  the  pressure  estimates  obtained 
at  both  measurement  surfaces  using  the  combination  projection  tech- 
niqe.  The  total  relative  error  in  the  pressure  estimations  (i.e.,  the  ab¬ 
solute  value  of  the  error  in  the  pressure  value  divided  by  the  pressure 
value)  that  were  obtained  for  the  2850  Hz  source  and  5000  Hz  source 
across  both  measurement  surfaces  are  shown  in  Table  1 .  The  com¬ 
bination  projection  technique  again  demonstrates  a  clear  advantage. 


Frequency 

Rigid  Array 

Open  Array 

Combination 

2800  Hz 

86.30 

184.01 

15.58 

5000  Hz 

36.16 

55.91 

30.43 

Table  1.  The  total  relative  error  in  free-field  pressure  estimation  is 
shown.  Details  as  in  text. 


4.  CONCLUSION 

In  this  paper,  we  present  a  novel  method  for  NAH  using  a  dual,  con¬ 
centric  SMA.  This  method  combines  holographic  projection  coeffi¬ 
cients  from  the  open  and  rigid  array  based  on  the  modal  coefficients. 
Field  measurements  demonstrate  that  this  technique  improves  the 
holographic  projection  of  the  acoustic  intensity  field.  In  future  work, 
we  will  examine  partial  field  decomposition  techniques  using  a  dual, 
concentric  SMA  (e.g.,  see  [12]). 
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Fig.  4.  Figure  4(a)  contrasts  combination  projection  with  outward 
naive  projection  from  the  inner  rigid  array,  and  Figure  4(b)  contrasts 
combination  projection  with  inward  naive  projection  from  the  outer 
open  array  (2850  Hz,  3.13  cm). 
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